DECOMPOSITION OF MONOMIAL ALGEBRAS: 
APPLICATIONS AND ALGORITHMS 

JANKO BOHM, DAVID EISENBUD, AND MAX J. NITSCHE 

Abstract. Considering finite extensions K[A] C K[B] of positive affine semigroup rings over 
a field K we have developed in [T] an algorithm to decompose K[B] as a direct sum of mono- 
mial ideals in if [A]. By computing the regularity of homogeneous semigroup rings from the 
decomposition we have confirmed the Eisenbud-Goto conjecture in a range of new cases not 
tractable by standard methods. Here we first illustrate this technique and its implementation 
in our Macaulay2 package MonomialAlgebras by computing the decomposition and the 
regularity step by step for an explicit example. We then focus on ring-theoretic properties of 
simplicial semigroup rings. From the characterizations given in [T] we develop and prove ex- 
plicit algorithms testing properties like Buchsbaum, Cohen-Macaulay, Gorenstein, normal, and 
seminormal, all of which imply the Eisenbud-Goto conjecture. All algorithms are implemented 
in our Macaulay2 package. 



1. Introduction 

Let B be a positive affine semigroup, that is, B is a finitely generated subsemigroup of N m 
for some m. Let K be a field and K[B] the affine semigroup ring associated to B, which can be 
identified with the subring of K [ti, . . . ,t m ] generated by monomials t u := t" 1 • . . . • i^ m , where 
u = (ui, . . . , u m ) G B. Denote by C(B) and by G(B) the cone and the group generated by B. 
From now on let A C B be positive affine semigroups with C(A) = C(B). We will now discuss 
the decomposition of K[B] into a direct sum of monomial ideals in -K"L4]. Observe that 

K[B] = Q) geG K.{t b \b£Bng}, 

where G := G{B)/G{A). Note that C(A) = C(B) if and only if K[B] is a finitely generated 
K [A] -module. From this it follows that G is finite, and we can compute the above decomposition 
since all summands are finitely generated. Moreover, there are shifts h g G G(B) such that 

I g :=K- {t h - h ° | b G B n 5 } 

is a monomial ideal in /fL4]. Thus, 

as Z m -graded K [A] -modules (with degt b = b). A detailed formulation of the algorithm comput- 
ing the ideals I g and shifts h g and a more general version of the decomposition in the setup of 
cancellative abelian semigroup rings over an integral domain can be found in [TJ Algorithm 1, 
Theorem 2.1]. 

Our original motivation for developing this decomposition was to provide a fast algorithm 
to compute the Castelnuovo-Mumford regularity reg K[B] of a homogeneous semigroup ring in 
order to test the Eisenbud-Goto conjecture [3]. Recall that the Castelnuovo-Mumford regularity 
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regM of a finitely generated graded module M over a standard graded polynomial ring R = 
X\, . . . ,x n ] is defined as the smallest integer m such that every j-th syzygy module of M is 
generated by elements of degree < m + j. Moreover, B is called a homogeneous semigroup if 
there exists a group homomorphism deg : G{B) — > Z with deg 6, = 1 for i = l,...,n, where 
Hilb(i?) = {61, . . . , b n } is the minimal generating set of B; by reg K[B] we mean its regularity 
with respect to the i?-module structure which is given by the -KT-algebra homomorphism R -» 
K[B], Xi 1 — ^ t b \ 

The toric Eisenbud-Goto conjecture can be formulated as follows: let K be a field and B a 
homogeneous semigroup, then reg K [B] < deg K[B] — codim K[B], where deg if [U] denotes the 
degree and codim K[B] := dimx K[B]± — dim K [B] the codimension. Even this special case of 
the Eisenbud-Goto conjecture is largely open, for references on known results see [TJ Section 
4]. The regularity of K [B] is usually computed from a minimal graded free resolution. If n 
is large this computation is very expensive, and hence it is impossible to test the conjecture 
systematically in high codimension using this method. However, choosing A to be generated by 
minimal generators e\, . . . , e<2 of C{B) of degree 1 the regularity can be computed as 

reg K [B] = max{regl 9 + degh g \ g G G}, 

where reg/ g denotes the regularity of I g with respect to the canonical T = K[x%, . . . , x^]-module 
structure given by T -» K [A], x^ 1— > t £i . Since the free resolution of every ideal I g appearing has 
length at most d—1, this computation is typically much faster than the traditional approaches. 
This enabled us to test the conjecture for a large class of homogeneous semigroup rings by using 
our regularity algorithm. See [TJ Section 4] for details. 

In Section [2j we illustrate, step by step, decomposition and regularity computation for an 
explicit example using our Macaulay2 [lj package MonomialAlgebras [2]. We say that K[B] 
is a simplicial semigroup ring if the cone C(B) is simplicial. In Section [3| we focus on simplicial 
semigroup rings K[B]. Based on the characterizations of ring-theoretic properties given in [lj 
Proposition 3.1] we develop explicit algorithms for testing whether K[B] is Buchsbaum, Cohen- 
Macaulay, Gorenstein, seminormal, or normal. We also discuss that, by known results, all these 
ring-theoretic properties imply the Eisenbud-Goto conjecture. The algorithms mentioned are 
implemented in our Macaulay2 package. 

2. Decomposition and regularity 

Our Macaulay2 package can be loaded by 
Macaulay2, version 1.4 

with packages: ConwayPolynomials , Elimination, IntegralClosure , LLLBases, 
PrimaryDecomposition, ReesAlgebra, TangentCone 
il : needsPackage "MonomialAlgebras"; 

We discuss the decomposition at the example of the homogeneous semigroup B C N 3 specified 
by a list of generators 

±2 : B = {{4, 0,0}, {2, 2,0}, {2, 0,2}, {0,2, 2}, {0,3,1}, {3, 1,0}, {1,1, 2}}; 
As an input for our algorithm we encode this data in a multigraded polynomial ring 

13 : K = ZZ/101; 

14 : S = K[x_l .. x_7, Degrees=>B] ; 
The command 

15 : dc = decomposeMonomialAlgebra S 

o5 = HashTable{ (-1,1,0) =>{ ideal (xi,x 3 ), (-1,1,0) } 

=>{ ideal 1, } } 

decomposes K[B] over K[A] where A C B is generated by minimal generators of C(B) with 
minimal coordinate sum; so in the example A = ((4, 0, 0), (2, 2, 0), (2, 0, 2), (0, 2, 2), (0, 3, 1)). The 
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keys of the hash table represent the elements of G and the values are the tuples (I g , h g ), hence 



(2.1) 



K[B] = (xT,x$ (-(-1,1,0)) ®K[A] 



as Z 3 -graded K[A]-modules; here we write K[A] = T/J with T = K[x\, X2, X3, x^, X5] and xl for 
the class of X{. Note that the on-screen output of MACAULAY2 does not distinguish between 
the class and the representative. To compute reg I g we will consider the standard grading on T: 



±6 
±7 
i8 
08 
i9 



o9 = 



o9 



KA = ring first first values dc; 

T = newRing(ring ideal KA,Degrees=>{5: l}) ; 

J = sub (ideal KA,T) ; 

ideal of T 



betti res 


total: 1 
1 
1 
2 

BettiTally 



The usual approach would be to obtain reg K[B] from a minimal graded free resolution of the 
toric ideal Ib with respect to the standard grading. 



ilO 
olO 
ill 
il2 



ol2 = 



ol2 



IB = monomialAlgebraldeal S; 
ideal of S 

R = newRing(ring IB ,Degrees=>{7 : l}) ; 
betti res sub (IB, R) 

1 2 3 4 5 
L 8 15 13 6 1 



total: 1 
1 
1 
2 
3 

BettiTally 



8 
3 
4 



10 6 1 



We observe that iegK[B] = 3. With degu = {u\ + 112 + us)/A by Equation (2.1) it holds 



reg K[B] = max {reg {x\,x%) + 0, mgK[A] + 0} 
By o9 we have reg-fT^A] = 2. To compute reg {x±, xs) do: 



il3 
ol3 
ol3 
il4 

ol5 
il5 



II = first (values dc)#0 
ideal (xi,x 3 ) 
ideal of KA 

g = matrix entries sub (gens II, T) ; 

1 2 
Matrix T < T 

betti res image map(coker gens J, source g, g) 
12 3 
2 



ol5 = 



5 
2 



4 1 



total : 
1 

2 .... 

3 .341 

Hence reg (xl, X3) = 3 and therefore reg K[B] = 3. Observe, that the resolution of K[B] has 
length 5, whereas the ideals I g have resolutions of length at most 3. The command 

il6 : regularityMA S 

0I6 : {3, {{ideal (xi,x 3 ) , (-1,1,0)}}} 
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provides an implementation of this approach, also returning the tuples (I g , h g ) where the maxi- 
mum is achieved. By [U Proposition 4.1] we have deg K[B] = j^G ■ degiT[A] = 10 since 

il7 : degree J 

ol7 = 5 

Moreover, codim K[B] = 4 since dimi-Tfi?] = dim C(B) = 3. Hence the ring K[B] satisfies the 
Eisenbud-Goto bound. 



3. Algorithms for ring theoretic properties 

In this section we focus on simplicial semigroup rings if[-B]. Based on the characterizations 
given in [TJ Proposition 3.1] we develop and prove explicit algorithms for testing whether K[B] is 
Buchsbaum, Cohen-Macaulay, Gorenstein, seminormal, or normal. Note that, in the simplicial 
case, all these properties are independent of K, and they imply the Eisenbud-Goto conjecture 
by results of El [6]. As an example, consider the following homogeneous simplicial semigroup 
BcN 3 specified by the generators 

il8 : B = {{4, 0,0}, {0,4,0}, {0,0, 4}, {1,0, 3}, {0,2, 2}, {3, 0,1}, {1,2,1}}; 
We compute the decomposition of K[B\ over K[A], where A = ((4,0,0), (0,4,0), (0,0,4)) C B 
is generated again by minimal generators of C(B) with minimal coordinate sum. 



il9 
i20 
o20 



S = K[x_l .. x_7, Degrees=>B]; 
decomposeMonomialAlgebra S 



HashTable{ 


(-1,0,1 


) =>{ ideal 1, 


(3,0 


1) 


} 




(-1,2,- 


1) =>{ ideal 1, 


(3,2 


3) 


} 




(0,2,2) 


=>{ ideal 1, 


(0,2 


2) 


} 




(1,0,-1 


) =>{ ideal 1, 


(1,0 


3) 


} 




(1,2,1) 


=>{ ideal 1, 


(1,2 


1) 


} 




(2,0,2) 


=>{ ideal (x 3 ,xi,x 2 ) , 


(2,0 


2) 


} 




(2,2,0) 


=>{ ideal 1, 


(2,2 


4) 


} 







=>{ ideal 1, 







} } 


K[B] 


^ k[A] e 


K[A](-1) A K[A](-2) 2 e 


(xi,x 2 , 


X3) 


(-1) 



Hence 



with respect to the standard grading induced by degu = {u\ + u 2 + «3)/4. It follows that 
depth K[B] = 1, thus, K[B] is not Cohen-Macaulay. Hence K[B] is also not normal by [5]. We 
can test seminormality via the following algorithm: 



Algorithm 1 Seminormality test 

Input: A simplicial semigroup B C N m . 
Output: true if K[B] is seminormal, false otherwise. 
1: Let ei, . . . ,e<z E B be minimal generators of C(B) with minimal coordinate sum, and set 
A := (ex, . . .,e d ). 

2: Compute B A ■= {x G B \ x £ B + (A \ {0})} as described in [1, Algorithm 1, Step 1]. 
3: for all x e Ba do 

4: Solve the linear system of equations Yli=i ^' e « = x for A = (Ai, . . . , Xd) G Q d - 
5: if IIAH^ > 1 then return false 
6: return true 



Here, by || — 1| we denote the maximum norm. Note that all Aj are non-negative since C(B) 
is a simplicial cone. Verifying in Step 5 the condition HA^ > 1 instead results in an algorithm 
which tests normality. Using our package we observe that K[B] is not seminormal: 

i21 : isSeminormalMA B 

o21 : false 
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The Buchsbaum property can be tested by the following algorithm. We denote by if [A] + the 
homogeneous maximal ideal of if [.A]. 



Algorithm 2 Buchsbaum test 

Input: A simplicial semigroup B = (b±, . . . ,b n ) C N m . 
Output: true if K[B] is Buchsbaum, false otherwise. 

1: Let ei,...,ed £ B be minimal generators of C(B) with minimal coordinate sum, and set 
A := (ei, . . .,e d ). 

2: Using the (minimal) generators e\, . . . , e d of A decompose 

where I g C K[A], h g G G(B) and G = G(B)/G{A) by [1, Algorithm 1]. 
3: if 3g £ G with I g / if [A] and I g / K[A] + then return false 
4: H := {h g \ g £ G with I g = K[A} + } 
5: C := {bi, . . . ,6„}\{0,ei,. . . ,e d } 
6: H + C := {h g + bi\h g e H,biG C} 
7: return true if (H + C) n H = and false otherwise. 



Proof. By jT] Proposition 3.1] the ring K[B] is Buchsbaum iff each ideal I g is either equal to 
if [A], or to if [A]+ and h g + 6 G B for all 6 G Hilb(i?). So, Step 3 is correct and we may now 
assume that I g = K[A] or I g = K[A} + for all g £ G. Recall that I g = {t v ~ h o \ v G T g }K[A], 
where F g = {x G -Ba I x G <?}. Moreover, note that | u G F g } is always a minimal 

generating set of I g and h g = Ylt=i mni {^fc I v e T 9 } e^, where v = Ylk=i K e k wrtn ^k G 

Since h + G -Ba for all /t G -ff and all A; = 1, . . . , d, we have H f] B = 0. In case that 
(if + C) n H ^ we obtain h + b ^ B for some h £ H and some b £ B \ {0}, that is, 
/i + Hiib(.B) ^ S. Hence if [B] is not Buchsbaum. 

In case that if [B] is not Buchsbaum, there is an h G H and some b G Hilb(-B) such that 
h + b ^ B. It is now sufficient to show that b G C and h + b G ii. By the above argument, 
6 G C. Let mfc = h + b + e& for k = 1, . . . , d. Suppose that m« ^ -Ba for some i G {1, . . . , d}. 
Since — et ^ B for all k = 1, ... ,d, necessarily rm — ej G B for some j ^ i. Consider 
y = rrij — Ylt=i n k^k £ B with G N such that 5^fc=i n fc * s maximal. By construction y G -Ba, 
moreover, rij = since — ej ^ -B. In the same way if x = rrii — ej — Ylt=i n k e k £ B with 
Ylt=i n k maximal, then x G Ba- Since rrii,mj G g for some g £ G, we also have x,y £ g. Since 
ei, . . . ,e<j are linearly independent, we have — A| > 2. Moreover, since t y ~ h9 , t x ~ hs £ K[A] 

we get that t y ~ hg is not a linear form. Hence I g ^ K[A] and I g ^ K[A] + , thus, £ Ba for all 
k = 1, . . . , d. We have #T g £ {1, d} by minimality, hence T g = {mi, . . . , m d }. By construction, 
h g = h + b and i 9 = if L4] + , therefore h + b £ H. □ 

Note that in Step 2 the shifts h g and hence the ideals I g are uniquely determined since 
ei, . . . ,ed are linearly independent. This is not true for arbitrary generating sets. By 
i22 : isBuchsbaumMA B 
o22 : true 

we conclude that K\B] satisfies the Eisenbud-Goto conjecture by [7]. Note that we can read 
off from the decomposition the regularity and the Eisenbud-Goto bound: we have regif [A] = 
and reg (xi, X2,x^) = 1, therefore regif [B] = max{0, 1, 2, 1 + 1} = 2. Moreover, degif [B] is the 
number of ideals which occur in the decomposition, hence degif [B] — codimif [B] =8 — 4 = 4. 

Note that, in case B is Buchsbaum, the regularity of K[B] is independent of the field if since 
all ideals in the decomposition are equal to the homogeneous maximal ideal or to if [A]. 

We finish this section by providing an algorithm for testing the Gorenstein property. 
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Algorithm 3 Gorenstein test 

Input: A simplicial semigroup B C N m . 

Output: true if K[B] is Gorenstein, false otherwise. 

1: Let ei, . . . , erf G B be minimal generators of C(B) with minimal coordinate sum, and set 
A := (ei, . . . ,e d ). 

2: Using the (minimal) generators e\, . . . , ej, of A decompose 

where I g C K[A], h g G G(B) and G = G(B)/G(A) by [1, Algorithm 1]. 
3: if 3g G G with J s ^ If [A] then return false 

4: := | g € G} 

5: if h & H with maximal coordinate sum is not unique then return false 
6: Let h € H with maximal coordinate sum. 
7: while H / do 

8: Let h g £ H 

9: if h — h g ^ H then return false 

10: # := #\ {/l,,, fc-ftfl} 
11: return true 



Proof. By jT] Proposition 3.1] the ring K[B] is Gorenstein iff I g = K[A] for all g G G and -ff 
has a unique maximal element with respect to < given by x < y if there is a z G 5 such that 
x + 2 = y. Note that H = Ba since I g = K[A] for all g G G. If there is a maximal element 
h G -ff , then this element has maximal coordinate sum. If H has more than one element with 
maximal coordinate sum, then H does not have a unique maximal element. To complete the 
proof we need to show that an element h g G H satisfies h g < h iff h — h g G H . But this follows 
from the fact that if x ^ Ba then x + y ^ -Byi for all x,y £ B. □ 

Note that performing Steps 1-3 of Algorithm [3] (and returning true afterwards) gives a test 
for the Cohen-Macaulay property. 
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